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Abstract 


Between  October  1987  and  August  1990,  two  arrays  of  inverted  echo  sounders  were  deployed  in  the 
Gulf  Stream  northeast  of  Cape  Hatteras  as  part  of  the  SYNoptic  Ocean  Prediction  Experiment. 
The  “Inlet  Array”  consisted  of  9  inverted  echo  soimders  (lES).  Centered  at  74°VV,  the  Inlet  Array 
was  designed  to  measure  key  parameters  that  describe  the  Gulf  Stream  path  variabibty  near  Cape 
Hatteras.  The  large  “Central  Array”  of  24  lESs  was  centered  on  the  current  near  ahmu 

400  km  downstream  of  the  Inlet  Array.  Spanning  nearly  300  km  in  both  the  cross-stream  and 
downstream  directions,  the  Central  Array  was  designed  to  monitor  the  thermocline  structure  of  the 
Gulf  Stream  in  the  region  of  large  meanders  and  frequent  ring  interactions. 

Using  objective  analysis,  we  have  mapped  the  Gulf  Stream  thermal  field  measured  by  the  lESs 
in  the  Inlet  Array.  In  this  report,  the  objective  analysis  technique  is  described  and  the  mapping 
parameters  are  documented.  Daily  maps  of  the  thermocline  depth  field  are  presented  for  the  period 
14  October  1987  through  31  August  1990. 
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1  Introduction 

Although  the  Gulf  Stream  has  been  studied  more  than  any  other  current  system,  many  fun¬ 
damental  questions  concerning  the  dynamics  and  energy  balances  governing  its  meandering  still 
remain.  To  address  these  questions,  a  multi-investigator,  multi-institutional  research  effort,  enti¬ 
tled  the  SYNoptic  Ocean  Prediction  (SYNOP)  experiment,  was  undertaken.  To  accomplish  the 
objectives,  the  observational  phase  of  SYNOP  was  designed  such  that  the  results  could  eventually 
be  integrated  with  numerical  models  of  the  Gulf  Stream. 

As  one  part  of  the  SYNOP  initiative,  Randy  Watts  of  the  University  of  Rhode  Island  and  John 
Bane  of  University  of  North  Carolina  conducted  a  joint  field  experiment  in  the  region  northeast 
of  Cape  Hatteras,  NC,  where  the  instantaneous  path  of  the  Gulf  Stream  is  convoluted  by  large 
time-varying  meanders.  The  study  area  extended  from  Cape  Hatteras  on  the  west  to  67°W  on 
the  east.  In  this  region,  the  Gulf  Stream  can  shift  laterally  through  an  envelope  which  exceeds  the 
instantaneous  width  of  the  current  itself.  The  Gulf  Stream  is  closely  associated  with  the  surrounding 
flow  fields  to  the  north  and  south  where  rings  frequently  form  and  interact  with  the  current. 

Between  October  1987  and  August  1990,  our  field  program  consisted  of  two  moored  arrays. 
The  “Inlet  Array”  was  composed  of  9  inverted  echo  sounders  (lES)  and  5  deep  current  meters. 
Centered  at  74°W,  the  Inlet  Array  was  designed  to  measure  key  parameters  that  describe  the 
path  variability  near  Cape  Hatteras  as  the  Gulf  Stream  leaves  the  continental  margin  and  flows 
into  deeper  water.  The  large  “Central  Array”  of  24  lESs  was  centered  on  the  current  near  AS^W, 
about  400  km  downstream  of  the  Inlet  Array.  The  Central  Array  also  consisted  of  thirteen  tall 
current  meter  moorings,  that  reached  from  the  bottom  into  the  high-velocity  Gulf  Stream  core. 
In  addition,  twelve  lESs  in  the  interior  of  the  array  were  outfitted  with  bottom  pressure  gauges. 
Spanning  nearly  300  km  in  both  the  cross-stream  and  downstream  directions,  the  Central  Array 
was  designed  to  monitor  the  Gulf  Stream  thermocline  structure  in  this  region  of  large-amplitude 
meanders  and  ring  interactions.  The  field  phase  of  the  program  ended  in  August  1990. 

The  subsequent  data  processing  tasks  were  divided  among  the  groups  at  URI  and  UNC.  The 
basic  steps  for  both  the  current  meters  and  lESs  included  transcription,  editing  and  conversion 
into  scientific  units.  The  lES  and  bottom  pressure  data  from  both  the  Inlet  and  Central  Arravs 
were  processed  at  URI.  These  data  have  been  documented  in  three  technical  reports  [Qian  et  ai, 
1990;  Fields  and  Watts,  1991a;  Fields  and  Watts,  1990b].  The  current  meter  records  from  the  Inlet 
Array  will  be  documented  in  Pickart  et  al.  [1991].  The  current  meter  records  for  the  Central  Array, 


processed  at  UNC,  will  be  documented  separately  [Shay  et  ai,  1991]. 

After  the  individual  lES  records  were  scaled  into  thermocline  depths,  we  combined  the  data  for 
each  array  to  produce  maps  of  the  Gulf  Stream  and  adjacent  thermal  fields.  This  report  documents 
the  mapping  procedures  for  the  Inlet  Array  for  October  1987  to  August  1990.  .Maps  of  the  Central 
Array  for  this  same  time  period  are  documented  separately  [Tracey  and  Watts,  1991], 
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2  Inlet  Array  lES  Data 

From  October  1987  to  August  1990,  the  Inlet  Array  consisted  of  9  lESs  located  on  three  cross- 
streatm  lines  (Figure  1).  During  the  summers  of  1988  and  1989,  all  instruments  were  recovered, 


longitude  (W) 

76  75  74  73  72  71  70  69  63  67 


Figure  1:  Inlet  and  Central  Arrays;  October  1987  to  August  1990.  IBS  sites  are  denoted  by  the 
open  circles.  The  box  outlines  the  140  km  by  120  km  region  of  the  Inlet  Array  which  has  been 
mapped  by  objective  analysis.  The  dashed  curve  indicates  the  mean  path  of  the  Gulf  Stream  for 
1975  to  1986  from  Gilman  and  Cornillon  [1990].  Bottom  depths  of  100,  200,  1000,  2000  m  are 
shown. 

refurbished,  and  redeployed.  The  IBS  positions  and  durations  for  the  three  deployment  periods 
are  given  in  Tables  1-3.  Unfortunately  during  the  second  deployment  period,  both  IBSs  ailong  the 
A  line  failed;  thus  the  array  was  limited  to  only  two  cross-stream  lines  from  June  1988  to  August 
1989. 

The  mean  downstream  distance  between  lines  A,  B,  and  C  was  55  km.  Along  the  A  and  C 
lines,  the  cross-stream  separation  distance  between  IBS  sites  was  approximately  40  km,  whereas 


on  the  central  B  line,  the  spacing  was  only  about  25  km.  The  array  wa.s  designed  so  that  key  path 
parameters,  displacement,  angle,  and  curvature,  could  be  measured.  Additionally,  the  lEb  loratlons 
were  selected  so  that  the  objectively  estimated  error  fields  would  indicate  accurate  mapping  of  the 
thermocline  topography  even  if  isolated  instruments  in  the  array  either  malfunctioned  or  were  lost. 
Furthermore,  a  couple  of  lES  sites  were  chosen  in  order  to  locate  them  along  the  suborbital  ground 
tracks  of  GEOSAT. 

After  recovery,  the  lES  data  were  edited  and  calibrated  using  a  suite  of  processing  steps  [Fields 
et  al,  1991].  Using  the  techniques  described  in  Watts  and  Johns  [1982],  Tracey  and  Watts  [1986] 
and  Howden  et  al.  [1991],  the  travel  times  measured  by  the  lESs  were  scaled  to  thermocline  depths 
and  subsequently  smoothed  using  a  40-huur  low-pass  filter  to  remove  the  inerticJ  and  tidal  signals. 
For  convenience,  the  thermocline  depth  is  expressed  as  the  depth,  in  meters,  of  the  12°C  isotherm 
(Z12).  The  data  were  subsampled  at  daily  intervals  prior  to  performing  the  objective  analysis. 
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Table  1:  IES88  Inlet  Array  -Locations  and  Duration 


lES  Site 

Latitude  (N) 

Longitude  (W) 

Start  Date 

End  Date 

IES88A1 

35*  18.56 

74*  36.94 

18  Oct  1987 

28  May  1988 

IES88A2 

35*  02.05 

74*  12.29 

14  Oct  1987 

12  Jun  1988 

IES88B1 

35*  45.06 

74*  27.97 

12  Oct  1987 

26  May  1988 

IES88B2 

35*  36.71 

74*  14.30 

14  Oct  1987 

28  May  1988 

PIES88B3 

35*  28.81 

74*  02.84 

13  Oct  1987 

28  May  1988 

PIES88B4 

35*  20.74 

73*  50.97 

13  Oct  1987 

28  May  1988 

IES88B5a 

35*  11.98 

73*  40.01 

13  Oct  1987 

18  Oct  1987 

IES88B5b 

35*  12.07 

73*  40.09 

21  Oct  1987 

14  Nov  1987 

IES88C1 

36*  04.54 

73*  56.98 

14  Oct  1987 

29  May  1988 

IES88C2 

35*  45.93 

73*  33.55 

14  Oct  1987 

30  May  1988 

Table  2:  IES80  Inlet  Array  -  Locations  and  Duration 


lES  Site 

Latitude  (N) 

Longitude  (W) 

Start  Date 

End  Date 

IES89A1 

35*  18.61 

74*  36.85 

Instrument  lost 

IES89A2 

34*  58.06 

74*  24.96 

No  data 

IES89B1 

35*  45.10 

74*  28.00 

29  May  1988 

9  Jun  1989 

IES89B2 

35*  36.92 

74*  14.20 

31  May  1988 

9  Aug  1989 

IES89B3 

35*  28.96 

74*  02.63 

30  May  1988 

9  Aug  1989 

35*  2C.70 

73*  5!  00 

.30  May  19«» 

10  Jun  1989 

IES89B4 

35*  20.80 

73*  50.50 

12  Jun  1989 

10  Aug  1989 

IES89B5 

35*  12.06 

73*  40.01 

30  May  1988 

10  Jun  1989 

IES89B5 

35*  12.04 

73*  39.95 

12  Jun  1989 

10  Aug  1989 

IES89C1 

36*  04.67 

73*  56.90 

1  Jun  1988 

41  Jun  1989 

IES89C2 

35*  46.20 

73*  32.90 

1  Jun  1988 

10  Aug  1989 

lable  3:  IES90  Inlet  Array  -  Locations  and  Duration 


lES  Site 

Latitude  (N) 

Longitude  (W) 

Start  Date 

End  Date 

IES90A1 

35*  12.32 

74*  43.91 

14  Aug  1989 

1  Sep  1990 

IES90A2 

34*  58.18 

74*  24.53 

12  Aug  1989 

1  Sep  1990 

IES90B1 

35*  45.13 

74*  27.90 

11  Jun  1990 

31  Aug  1990 

1IES90B2 

35*  37.01 

74*  13.82 

Instrument  lost 

TIES90B3 

35*  30.07 

74*  03.40 

11  Aug  1990 

1  Sep  1990 

TIES90B4 

35*  20.75 

73*  50.60 

12  Aug  1990 

1  Sep  1990 

TIES90B5 

35*  12.13 

73*  39.66 

12  Aug  1990 

1  Sep  1990 

IES90C1 

36*  04.57 

73*  56.84 

13  Jun  1990 

2  Sep  1990 

TIES89C2 

35*  46.15 

73*  33.00 

13  Aug  1989 

16  Oct  1989 

TIES90C2 

35*  46.22 

73*  32.75 

18  Oct  1989 

2  Sep  1990 
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3  Objective  Mapping  of  the  Thermocline  Depth  Field 

In  order  to  produce  gridded  maps  of  the  Gulf  Stream  Z\2  field,  we  used  objective  analysis  (OA)  to 
interpolate  between  the  irregularly-spaced  observations.  The  interpolation  is  made  at  each  output 
grid  point  as  a  weighted  average  of  input  values  at  the  neighooring  measurement  sites.  The  optimal 
weighting  was  determined  according  to  the  Gauss-Markov  theorem  using  knowledge  of  the  space- 
time  correlation  function  (see  below).  For  our  application,  the  correlation  function  was  dependent 
on  radial  distance  and  time,  p(r,t'). 

The  basic  OA  equations  are  clearly  formulated  in  Bretherton  et  al.  [1976]  and  e.xtensions  to 
anisotropic  fields  are  given  in  Carter  [1983]  and  Carter  and  Robinson  [1987].  The  method  was 
further  modified  by  Watts  and  Tracey  [1985]  and  Watts  et  al.  [1989]  to  apply  the  technique  to  the 
Gulf  Stream  frontal  region.  To  summarize  briefly,  their  adaptations  include  preconditioning  the 
observations  by  removing  a  spatially-dependent  mean  thermal  field  and  subsequently  normalizing 
the  residual  perturbations  by  the  standard  deviation  (STD)  field.  These  techniques  combine  to 
produce  perturbation  fields  which  have  homogeneous  statistics,  a  requirement  of  objective  analysis. 
After  running  the  OA.  both  the  mean  and  STD  fields  are  restored  to  the  output  fields,  producing 
maps  of  the  Gulf  Stream  Zu  field.  We  used  the  procedures  of  Watts  et  al.  [1989]  to  generate  the 
OA  maps  shown  in  this  report. 

Since  we  have  a  large  set  of  observations,  we  adapted  programs  generously  provided  by  E.  Carter 
(documented  in  Carter  and  Robinson  [1987])  that  make  the  computing  much  more  efficient.  We 
took  advantage  of  the  facf  that  data  values  which  are  distant  (both  spatijiUy  and  temporally)  from 
tho  output  point  have  little  influence  on  the  estimate  when  compared  with  the  nearby  observations. 
We  increased  the  efficiency  by  restricting  the  available  data  to  a  smaller  subset  prior  to  running 
the  OA  at  each  output  point.  First,  we  specified  a  maximum  time  lag  (tJnax)  maximum  radial 
separation  distance  (rm„)  and  eliminated  observations  which  exceeded  these  restrictions.  For  the 
Inlet  Array  OA  maps,  we  specified  =  ±1  day  and  =  120  km.  When  determining  the 
radial  distance  for  each  data  point,  we  took  into  account  downstream  phase  propagation  of  Gulf 
Stream  meanders  (discussed  below)  such  that 

r=  [(z'-cO'  +  (yy]'  (1) 

where  t.ie  phase  speed  c  =  ?4  km  d~*. 

To  further  reduce  the  data  around  a  given  output  grid  point  to  use  only  the  most  influential 
input  data,  we  selected  a  smaller  subset  of  these  restricted  observations.  We  chose  only  the  .V  points 


witli  the  highest  correlations  to  be  used  for  the  estimation  at  the  output  grid  point.  The  choir.' 
for  iV  is  a  trade  off  between  wanting  enough  observations  to  reduce  random  noise  and  wanting 
fewer  points  to  increase  computational  efficiency.  The  Gauss-Markov  procedure  involves  inversion 
of  an  .V  X  iV  matri.x.  repeated  for  each  output  grid  point  of  every  map.  If  N  is  too  large  (such 
that  the  input  data  are  themselves  not  sufficiently  independent)  then  the  m.atrix  inversions  arr 
computationally  unstable  and  the  errors  actually  increase.  We  used  \  -  6  points  to  produce  the 
Inlet  Array  Z12  maps. 

While  the  choices  for  ^m»x.  ^nd  \  place  restrictions  on  the  data  to  be  used  for  inter¬ 

polation,  they  do  not  directly  appear  in  the  OA  equations.  Aside  from  the  correlation  function, 
the  only  user-specified  parameter  that  is  explicitly  a  factor  in  the  equations  is  the  noise  level  of 
the  data  (c).  By  specifying  a  non-zero  noise  level,  the  interpolation  scheme  is  not  required  to  fit 
the  data  values  exactly,  and  some  smoothing  is  permitted.  We  specified  c  =  0.05  (in  terms  of 
percentage  of  the  data  variance),  since  it  is  our  best  estimate  of  the  true  error  associated  with  the 
lES  measurements. 

3.1  The  Space-Time  Correlation  Function 

To  produce  the  OA  maps  of  the  Inlet  Array  thermal  field,  it  wa-  necessary  to  determine  a  correlation 
function  appropriate  for  that  region  of  the  Gulf  Stream.  We  followed  the  procedures  described  in 
Watts  et  al.  [1989]  to  calculate  the  spatial  and  temporal  correlations  from  our  observations.  In 
all,  we  used  Z12  records  from  28  lESs  obtained  during  three  separate  deployment  periods  between 
1987  and  1990.  We  csdculated  the  spaticd  separations  {x'.y')  between  pairs  of  lESs  within  each 
deployment  period  and  determined  their  cross-correlations  for  several  time  lags  t' .  Subsequently, 
the  correlations  were  smoothed  by  a  60-km  Gaussian-weighted  low-pass  filter  to  remove  some  of 
the  small  scale  noise.  Figure  2  shows  the  observed  correlations  in  the  full  {x',  y')  plane  for  lag  times 
of  0,  1,  and  2  days. 

Next,  to  examine  the  observed  correlations  fc  spatial  anisotropy,  we  sorted  the  pairs  into 
downstream,  cross-stream  and  diagonal  categories  and  recomputed  the  correlations  for  t'  =  0  days. 
Figure  3  reveals  the  similarity  between  the  three  groups,  indicating  no  significant  directional  de¬ 
pendence.  Thus,  we  expressed  the  spatial  decay  of  the  correlations  as  a  function  of  radial  distance, 
rather  than  of  both  x'  and  y' .  The  best  fit  to  the  observed  correlations,  shown  in  Figure  3,  was  an 
exponential-cosine  function  of  the  form 

p{r)  =  Foexp  fos 


(right)  ol 


where  =  (z')^  +  iv'Y,  A  =  254  km,  and  B  =  160  km.  The  scale  factor  is  defined  as  fo  =  1  - 
where  c  =  0.05  is  the  error  assoicated  with  the  observations. 

To  determine  the  temporal  decay  of  the  correlation  function,  we  calculated  the  auto-correlations 
of  19  lESs  using  positive  time  lags  of  0  to  24  days.  For  this  analysis,  we  were  careful  to  select 
only  lESs  with  long  deployment  periods,  requiring  the  record  length  to  exceed  180  days.  The 
autocorrelations,  averaged  for  each  time  lag,  are  shown  in  Figure  4.  To  determine  an  analytic 
expression  for  the  temporal  decay  rate,  we  fitted  an  exponential-cosine  of  the  form 

p(t')  =  exp(-ait')cos  (3) 

to  the  correlations.  From  symmetry  requirements,  the  correlations  for  negative  time  lags  are 
p(-x\—y',—i')  =  For  the  full  set  of  correlations  (time  lags  out  to  24  days),  the 

best  fit  is  shown  by  the  dashed  curve  in  Figure  4  and  is  obtained  by  specifying  ai  =  0.0527  d“'  and 
Oj  =  —2.3056  d“*.  However,  this  curve  does  not  adequately  describe  the  rapid  decay  rate  of  the 
shorter  time  lags.  Since  we  only  use  lags  out  to  2  days  when  performing  the  objective  mapping,  we 
wanted  to  use  a  temporal  decay  rate  which  is  more  appropriate  for  these  shorter  time  scales.  By 
specifying  aj  =  0.1182  d”*  and  aj  =  0.0831  d“*,  we  obtained  the  solid  curve  shown  in  Figure  4. 
This  is  the  temporal  decay  rate  used  to  define  the  correlation  function  for  the  Inlet  Array. 

Subsequently,  we  examined  the  propagation  of  the  correlation  peak  with  time.  We  expected  the 
peak  to  move  downstream  at  a  rate  comparable  to  that  of  Gulf  Stream  meanders.  We  selected  21 
p2urs  of  lESs  which  were  located  downstream  of  one  another  and  cjilculated  their  cross-correlations 
for  t'  =  0  to  20  days.  For  each  pair,  we  found  the  time  associated  with  the  maximum  correlation 
and  determined  the  speed  based  on  the  spatial  separation  between  the  instruments.  The  rates 
ranged  from  10-55  km  d“',  averaging  43  kra  d”*  ±  14  km  d~’ .  However,  in  Figure  2,  the  observed 
correlation  peak  at  1'  =  2  days  occurs  at  x'  =40-60  km,  indicating  phase  speeds  of  only  20- 
30  km  d~^.  To  examine  the  effect  of  phase  speed  (c)  on  the  output  Z12  maps,  we  produced  three 
sets  of  maps  using  c  =  12,  24,  and  42  kra  d"’.  Then  we  determined  the  mean  and  standard 
deviation  of  the  differences  between  the  mapped  Zu  fields.  For  the  18  Zjj  maps  examined,  there 
were  only  small  differences  between  the  maps  generated  with  c  =  12  km  d"!^  and  c  =  24  km  d~*. 
The  mesm  difference  was  less  than  one  meter  and  the  standard  deviation  was  under  3  m.  Not 
surprisingly,  the  largest  differences  occurred  between  the  maps  produced  with  c  =  12  km  d~’  and 
c  =  42  km  d“* .  But  even  then,  the  mean  difference  between  the  18  OA  maps  was  less  than  4  m  and 
the  standard  deviation  was  about  7  m.  In  all  cases,  the  largest  differences  (about  30  m)  occurred  in 


Distance  (km) 

Figure  3:  Correlation  function  for  downstream,  cross-stream,  and  diagonal  lES  pairs  versus  radial 
distance.  The  solid  curve  is  the  analytic  function  p{r). 


5.0  10  0  15-0  20. C 

Time  Lag  (days) 


Figure  4;  Average  auto-correlations  for  19  lE.Ss  versus  lag  time.  The  solid  and  da.sliod  (uim-s 
are  the  best  fits  of  the  analytic  function  p[t')  to  the  correlations  at  shorter  and  longer  time  la-s, 
respectively. 
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the  extrapolated  portions  of  the  mapping  region.  Since  the  average  differences  between  the  maps 
are  quite  small  compared  with  the  mapping  accuracy  of  47  m,  the  choice  for  phase  speed  is  not 
crucial  to  the  quality  of  the  output  maps.  We  chose  a  mid-range  phase  speed  of  c  =  24  km  d“'  (cf. 
Equation  [1])  to  produce  the  Inlet  Array  OA  maps. 

Finally  we  combined  Equations  [2]  and  [3]  to  obtain  the  full  analytic  expression  for  the  corre¬ 
lation  function  of  the  thermocline  depth  in  the  Inlet  Array: 

=  F’oexp(-oit')cos  exp  cos  (4) 

where  r  is  defined  in  Equation  [1].  Figure  2  shows  piniet  in  plan  view  for  the  same  time  lags  as  the 
observed  correlations. 

3.2  Mean  and  STD  Fields 

Watts  et  al.  [1989]  showed  that  the  quality  of  the  OA  maps  is  influenced  most  strongly  by  the  choice 
of  mean  Z12  field  that  is  removed  from  the  observations  prior  to  performing  the  objective  analysis 
(and  restored  later).  They  found  that  the  best  accuracy  was  obtained  with  a  time-averaged  mean 
field.  In  the  past  we  determined  the  temporal  mean  field  by  fitting  a  third-order  polynomial  to  the 
deployment-long  average  Zn  values  obtained  at  the  lES  sites.  Unfortunately,  this  method  tended 
to  give  unrealistic  values  along  the  edges  of  the  mapping  region  and  we  were  required  to  tweek  the 
inputs  to  obtain  a  smooth  field.  Thus  for  the  Inlet  Array,  we  elected  to  determine  the  mean  Z12 
field  manually.  First,  we  combined  the  Z12  averages  for  all  lES  sites  from  the  three  deployment 
periods  onto  a  single  map.  Next,  we  picked,  by  eye,  the  relative  cross-stream  positions  of  several 
thermocline  depths  along  the  upstream  and  downstream  edges  of  the  mapping  region.  Then  we 
used  a  cubic  spline  to  interpolate  between  these  positions  to  obtain  Z12  values  at  20  km  intervals 
along  the  edges.  Finally,  we  linearly  interpolated  between  the  two  edges  to  fill  in  the  remainder  of 
the  grid.  The  mean  Zu  field  chosen  is  shown  in  Figure  5. 

The  standard  deviations  of  the  28  Inlet  Array  lESs  varied  systematically  with  the  mean  Z12 
value,  related  to  the  proximity  of  the  lES  sites  to  the  mean  location  of  the  steeply  sloping  ther¬ 
mocline.  This  relationship  is  shown  in  Figure  6.  We  approximated  the  STD  field  by  a  Gaussian 
function  of  the  form 

Inlet  =  D  +  Cc~^ 

where  D  =  16  m  and  C  =  80  m.  The  exponent  is  defined  as  T  =  ^  where  A  =  470  m 

and  B  =  270  m.  The  STD  field  for  the  Inlet  Array  is  shown  in  Figure  5  along  with  the  associated 


-<0  0  *0  80  120  -40  0  40  80  120 

Distance  (km)  Distance  (km) 

Mean  Standard  Deviation 

Figure  5:  Mean  Zu  (left)  and  standard  deviation  (right)  fields  are  contoured  in  plan  view.  Contour 
intervals  are  50  m  for  the  mean  field  and  10  m  for  the  standard  deviation  field.  lES  sites  are  indicated 
by  the  +  marks.  Axis  tick  marks  are  at  20  km  intervals. 

mean  Zu  field. 

3.3  Construction  of  the  Zu  Maps 

The  objective  analysis  was  performed  on  the  “perturbation”  lES  Zn  records  (i.e.,  the  means  were 
removed  and  the  standard  deviations  normalized).  The  correlation  function  piniet  was  computed 
using  Equation  [4]  each  time  a  {x',y',t')  triplet  was  specified.  This  resulted  in  a  smooth,  continu¬ 
ously  changing  function  which  produced  better  quality  output  maps  than  the  lookup  table  method 
that  we  had  used  previously.  Data  restrictions  were  controlled  by  several  user-specified  parameters, 
summarized  in  Table  4. 

Table  4:  User-specified  control  parameters  for  OA  mapping. 


Parmeter 

Specification 

N 

6  points 

^mar 

i20  km 

f'mox 

±1  day 

C 

0.05 

c 

24  km  d”’ 

Observed  STD  (in 


i.i 


Figure  6:  Observed  standard  deviations  for  28  lESs  during  three  deployment  periods  from  1987- 
1990  versus  mean  Zu  depth.  The  solid  curve  is  the  analytic  function  crioiet- 

The  output  of  the  objective  mapping  is  the  perturbation  field  on  a  full  grid  of  points,  with 
20  km  spsicing.  The  perturbations  were  converted  to  thermocline  depths  by  renormalizing  by  the 
STD  field  and  restoring  the  mean  Z12  field.  Subsequently  the  Z\2  fields  were  smoothed  using  a 
second-order  Shapiro  filter  to  remove  small-scale  (~  2Ai  waves)  features  and  the  outermost  grid 
points  were  dropped.  The  final  size  of  the  Z12  grid  for  the  Inlet  Array  was  140  km  by  120  km  (8 
by  7  grid  points).  From  June  1988  to  August  1989,  when  no  data  were  obtained  along  the  A  line, 
the  mapping  region  was  shortened  to  80  km  by  120  km  (5  by  7  points). 


4  Estimated  Error  Fields 


Statistical  estimates  (percent  standard  deviation)  of  the  accuracy  of  the  output  fields  were 
produced  during  the  objective  analysis.  The  error  associated  with  each  output  grid  point  depends 
only  on  the  locations,  in  space  and  time,  of  the  observations  and  the  correlation  function.  The  error 
is  independent  of  the  measurements  themselves.  Thus,  the  same  error  field  applies  for  all  daily  OA 
maps  produced  from  a  given  array  of  instruments. 


Figure  7:  Estimated  error  (percent  STD)  fields  for  four  different  array  configurations  apply  to  the 
Z\2  fields  for  the  dates  shown.  Errors  greater  than  0.18  are  contoured  at  intervals  of  0.03.  For 
clarity,  errors  between  0.30  and  0.48  are  indicated  by  datshed  contours.  lES  sites  are  shown  by  the 
■  +  '  marks.  The  common  plotting  grid  was  chosen  for  the  four  arrays  out  of  convenience. 

During  the  34-month-long  deployment  period,  there  were  four  different  array  configurations  due 
to  the  loss  or  repositioning  of  several  lESs.  The  corresponding  error  fields,  along  with  the  lES  sites, 
arp  shown  in  Figure  7.  The  four  error  fields  are  basically  the  same  because  the  close  spacing  of  the 
1F.S^  gavr  rise  to  low  errors  throughout  the  mapping  region  even  when  isolated  instruments  were 


lost.  The  errors  are  low  near  the  lES  sites  and  increase  with  distance  from  those  sites.  During 
the  cruises  to  cycle  the  array  (May-June  1989,  May-June  1990,  August  1989,  and  .4ugust  1990), 
the  array  configurations  changed  frequently  while  the  lESs  were  being  recovered  and  redeployed. 
During  these  periods,  the  actual  error  fields  associated  with  the  Z\i  maps  changed  daily.  These 
‘transitional'  error  fields  are  not  shown  in  Figure  7,  but  areas  of  higher  error  are  indicated  on  the 
individual  objective  maps. 

We  use  the  estimated  error  fields  to  mask  out  portions  of  the  OA  maps  where  the  map  quality 
is  predicted  to  be  poor.  On  the  Z\i  maps,  regions  with  predicted  error  30-48%  are  shaded  by 
dashed  contours,  and  larger  errors  are  indicated  by  solid  contours.  Based  upon  analogies  drawn 
from  Watts  et  al.  [1989],  our  preliminary  estimate  is  that  the  dimensional  uncertainty  of  the  Z\2 
field  is  less  than  50  m  within  the  mapping  region  for  which  the  predicted  error  is  less  than  30%. 
We  have  used  50  m  as  the  contouring  -nterval  for  the  subsequent  maps  of  the  Z12  field. 


5  Average  Z12  Maps 


To  produce  average  maps  of  the  Z12  field  in  the  Inlet  Array,  we  summed  the  maps  (presented  in  the 
following  section)  by  grid  point  and  then  divided  by  the  number  of  maps.  We  averaged  two  ways, 
by  deployment  period  and  by  season.  The  deployment-long  average  Z12  fields  for  the  three  periods 
are  shown  in  Figure  8  with  the  averaging  dates  indicated  below  each  map.  Next,  we  divided  the 
maps  into  twelve  groups  of  approximately  90  days,  corresponding  to  the  four  seasons.  Two  of  the 
groups  were  for  shorter  time  periods  due  to  the  timing  of  the  start  and  end  dates;  the  group  for 
October- December  1987  consisted  of  81  days  and  the  July-September  1990  group  of  only  63  days. 
Average  Z12  maps  were  produced  for  each  group  and  these  are  shown  in  Figure  9. 
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Figure  8:  Z12  fields  averaged,  by  grid  points,  for  the  three  deployment  periods.  Averaging  peri¬ 
ods  are  indicated  below  each  map.  Depth  contours  are  in  meters.  Distances  on  the  axes  are  in 
kilometers. 
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6  Daily  Z\2  Maps  for  the  Inlet  Array 

Objective  maps  of  the  Gulf  Stream  thermal  field  were  produced  at  daily  intervals  for  October  14, 
1987  through  August  31,  1990.  The  maps  are  shown  for  1200  UT  on  the  dates  indicated  below 
each  frame. 

The  frames  correspond  to  the  boxed  region  shown  in  Figure  1.  The  actual  lES  sites  are  indicated 
by  the  ‘+’  marks  and  the  positions  are  listed  in  Tables  1-3.  Each  lES  position  was  converted  to  grid 
coordinates  using  scale  factors  of  111  km  per  degree  latitude  and  86.263  km  per  degree  longitude; 
distances  were  referenced  to  the  grid  origin  at  35°N,  74° W.  Subsequently  the  grid  was  rotated 
50°  so  that  the  base  was  oriented  040°T.  This  rotation  was  done  prior  to  running  the  objective 
analysis  for  plotting  convenience,  not  as  a  requirement  of  the  OA  technique.  The  latitude-longitude 
positions  of  the  grid  corners  are  listed  in  Table  5. 

The  Z\2  field  is  contoured  at  50  m  intervals.  Extrapolated  regions  of  the  OA  maps  have  been 
masked  out  where  the  the  estimated  errors  are  high;  errors  of  30-48%  are  indicated  by  dashed 
contours  and  higher  errors  by  solid  contours.  The  contour  interval  of  the  estimated  errors  is  3%. 


Table  5:  Corner  Positions  of  Mapping  Regions 


Map  Dates: 

14  October  1987  to  12  June  1988 
and 

10  August  1989  to  31  August  1990 

Corner 

X  (km) 

y  (km) 

Latitude  (N) 

Longitude  (W) 

Upper  Left 

-40 

+  100 

35°  18.18 

75°  11.17 

Lower  Left 

-40 

-20 

34°  36.49 

74°  07.23 

Upper  Right 

-1-100 

+  100 

36°  16.15 

74°  08.57 

Lower  Right 

-1-100 

-20 

35°  34.46 

73°  04.63 

Map  Dates: 

13  June  1988  to  9  August  1989 

Corner 

I  (km) 

y  (km) 

Latitude  (N) 

Longitude  (W) 

Upper  Left 

-1-20 

100 

35°  43.03 

74°  44.34 

Lower  Left 

-1-20 

-20 

35°  01.33 

73°  40.40 

Upper  Right 

+  100 

+  100 

36°  16.15 

74°  08.57 

Lower  Right 

+  100 

-20 

35°  34.46 
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